Rydberg trimers and excited dimers bound by internal quantum reflection 
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Quantum reflection is a pure wave phenomena that predicts reflection of a particle at a changing 
potential for cases where complete transmission occurs classically. For a chemical bond, we find that 
this effect can lead to non-classical vibrational turning points and bound states at extremely large 
interatomic distances. Only recently has the existence of such ultralong-range Rydberg molecules 
been demonstrated experimentally. Here, we identify a broad range of molecular lines, most of 
which are shown to originate from two different novel sources: a single-photon associated triatomic 
molecule formed by a Rydberg atom and two ground state atoms and a series of excited dimer states 
that are bound by a so far unexplored mechanism based on internal quantum reflection at a steep 
potential drop. The properties of the Rydberg molecules identified in this work qualify them as 
prototypes for a new type of chemistry at ultracold temperatures. 



The development of techniques to manipulate ultra- 
cold atoms has opened up new avenues for molecular 
physics. The recent experimental success in producing 
cold ground state molecules [TJ [2j [3] is paving the way 
towards studying chemistry at ultralow temperatures [4] . 
This regime also offers unprecedented opportunities for 
the control of large-scale molecules involving highly ex- 
cited states [5j [6] , which hold great promise for such ap- 
plications due to the high sensitivity of Rydberg atoms 
to external fields [7 . 

A particular class of Rydberg molecules arises from 
binding between a Rydberg and a ground state atom, 
which is based soley on low-energy collisions between the 
Rydberg electron and the ground state atom [8]. Com- 
pared to the large de Broglie wavelength of the electron, 
the ground state atom represents a point-like perturba- 
tion which probes directly the Rydberg atom's wavefunc- 
tion. Consequently, the molecular potential shows the os- 
cillatory character of the corresponding probability den- 
sity. Owing to the large size of Rydberg atoms, binding 
occurs at extremely long range on the order of several 
thousand Bohr radii, ao- 

Recently, we have provided experimental proof for the 
existence of these ultralong-range molecules in an ultra- 
cold rubidium gas [6]. A simple model with an adjustable 
s-wave scattering length was found to reproduce the bind- 
ing energies of the vibrational ground states for a range 
of principal quantum numbers. 

In this work, we demonstrate the existence of several 
additional molecular lines. The observed spectra can 
be adequately described with Fermi's pseudopotentials 
for local Rydberg electron scattering off the perturbing 
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ground-state atom [10] , but within a non-perturbative 
Green's function approach which accounts for the strong 
collisional couplings between Rydberg state manifolds. 
Despite the lack of an inner potential barrier, the calcu- 
lated excited states have the peculiar property to avoid 
smaller interatomic distances. As we will show, this is 
due to internal quantum reflection at a shape resonance 
of electron-atom scattering. In addition to reproducing 
all previously observed lines [6 , the calculations predict 
additional dimer and trimer states that are observed in 
our new experiments. 



I. THE MOLECULAR POTENTIAL 

The state of the molecular Rydberg electron is deter- 
mined by both long-range attraction to its parent ion 
as well as the short-range interaction with the nearby 
ground state atom. This vast length scale disparity 
poses considerable problems for current quantum chem- 
istry approaches, which can be elegantly resolved through 
Fermi's pseudopotential approximation [8 . We first ap- 
ply the Born-Oppenheimer approximation and a single- 
active particle treatment of the excited electron to obtain 
the following Hamiltonian (in atomic units) 

H = + V eA+ (|r - R|) + V eA (r) = A + V eA (r) (1) 

of the molecular Rydberg electron. Here, H A is the 
Hamiltonian of the unperturbed Rydberg atom, whose 
center is located at R, and where the molecular axis R 
also defines the quantization axis. V eA (r) accounts for the 
interaction of the Rydberg electron with the perturbing 
ground state atom, which is described by a Fermi-type 
pseudopotential for electron atom scattering [9j [10] , in- 
cluding the first two / = 0, 1 partial waves. 

Within the Fermi model [9], the molecular interaction 
of \ow-£ Rydberg states is commonly described by the 
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FIG. 1: Molecular potential curve for n — 35. Results ob- 
tained from the Green's function calculation (solid line) and 
from the model potential Eq.Q (dotted line) with an effective 



scattering length A, 



(eff) 



-19.48 clq. The grey area around 



R « 1200 do marks the innermost "badland" region defined by 
dX/dR > 1 [24]. The baseline of the molecular wave functions 
corresponds to their respective vibrational energies. The in- 
set shows the rapid potential drop at the avoided crossing and 
illustrates the variable inner boundary condition at a distance 
Ro. 



simple first-order perturbative potential [9j [10] 

UW(R) = 2nA s \<p n (R)\ 2 + 6nA 3 p \ W„(#)| 2 (2) 

where i(j n is the unperturbed electronic Rydberg wave 
function and A s ^(k) is the energy dependent s(p)-wave 
scattering length. The electron momentum k depends on 
the Rydberg energy E n through the quasiclassical rela- 
tion E n = k 2 /2 — 1/R. This approach yields an intrigu- 
ingly simple and intuitive description of the molecular in- 
teraction and is typically employed to describe molecules 
formed with energetically well isolated ns Rydberg states 
[HI [12]. Even for this case, however, Eq.Q turns out to 
be inaccurate at the level of accuracy of the experimen- 
tal data. Hence, we performed non-perturbative Green's 
function calculations [11] of the molecular potentials (see 
Methods). 

As an example, we compare the calculated potential 
energy curves U (R) with the first-order perturbative po- 
tential UM(R) for n = 35 in Fig.[lJ As is known for high- 
£ Rydberg molecules, Eq.Q fails in the vicinity of the 
p-wave shape resonance [13] [14] because the resonance 
leads to an avoided potential crossing. This curve cross- 
ing is crucial for determining the vibrational spectrum 
of the molecule found in our experiments. Unexpectedly, 
also falls short of more accurate calculations around 
the outer potential well at R = 1900 clq. 



One can, nevertheless, obtain an approximate descrip- 
tion of the strong molecular resonance, localized in the 
outermost minimum of the potential energy curve, by 
using an effective, properly adjusted s-wave scattering 
length A { s e ^\k = 0) with Eq.Q (see dotted line in 
Fig. [T]). The corresponding vibrational energy for all 
principal quantum numbers n = 35,36,37 is well re- 
produced for an effective scattering length of Ai e ^ = 
— 19.48 do. While this value is 20% larger than the actual 
value of A s (k = 0) (see Sec. II), it may serve as an effec- 
tive parameter or characteristic of a given atomic species 
that yields the respective n-dependence of the strongest 
molecular resonance, based on the simple model potential 



II. EXPERIMENTAL RESULTS 

To create the ultralong-range Rydberg molecules, we 
prepare a cold sample of 87 Rb atoms in the ground state 
5si/2j-^ = 2,77ip = 2. The nsi/ 2 Rydberg state of prin- 
cipal quantum numbers, n = 35 — 37, is addressed by 
two narrow-band lasers (see Methods). The vibrational 
bound dimer and trimer states can be created directly via 
photoassociation if the detuning of the excitation light 
from the atomic Rydberg state matches the molecular 
binding energy, Eb- Therefore, the molecular states ap- 
pear in the spectra on the red side of the atomic Ryd- 
berg line, as shown in Fig. [2] The uppermost graph (a) 
shows the full Rb(35s) spectrum which consists of two 
atomic states m s — +1/2 (j) and m s = —1/2 (j), sepa- 
rated by the Zeeman splitting A# in the magnetic trap, 
and several molecular lines. Due to the spin polariza- 
tion in the magnetic trap (m F = +2), only the molec- 
ular triplet states (||) are excited. The atomic state 
ns 1 / 2 ,m s = +1/2 and the molecular state 3 £(ns — 5s) 
have the same Zeeman shift (see inset) and thus the bind- 
ing energy Eb can be directly measured as the energy 
difference between the corresponding lines in the spec- 
trum. Note that in our previous work [6 the binding 
energies Eb were defined differently. Fig.p] (b)-(d) show 
the 35s to 37s spectra at high resolution. 

After field ionization, we observe both an atomic ion 
signal corresponding to Rb + and a molecular ion signal 
corresponding to RbJ which can be clearly identified in 
the time-of-flight spectra. The Rb + and RbJ spectra 
predominantly show the same features but with different 
intensity distributions. A few deep bound states exclu- 
sively appear in the RbJ spectrum, but can be assigned 
as dimer states by our theoretical model. The fact that 
the strengths of the two ion signals are comparable for 
the molecular states is a clear indication that RbJ is the 
product of an additional decay that occurs for the molec- 
ular states. To obtain maximum spectral information, we 
use both the Rb + as well as RbJ for the following anal- 
ysis of the molecular data. 

To assign the observed lines based on the described 
potential calculations, we employ the ab-initio data of 
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FIG. 2: Spectra for principal quantum numbers n — 35 — 37. 
a) Full spectrum for Rb(35s) and illustration of the atomic 
Zeeman shift A#. b) - d) Magnified spectra of the states 35s 
to 37s. Frequencies in all spectra are measured with respect 
to the atomic transition 5si/ 2 — > nsi/ 2 . The upper parts show 
the atomic ion spectrum Rb + (black) and the lower ones the 
spectra of the molecular ions Rb J (blue) . Calculated binding 
energies (Tab. [TJ> are indicated by vertical lines. 



|T2l [13] for the s-wave scattering phaseshift So(k). Since 
the available p-wave scattering data [I2j [15] does not 



yield an accurate description of the measured spectra, 
we use the modified effective range expression [16] for 
the corresponding phaseshift o\. This three-parameter 
fitting procedure allows the assignment of most of the 
observed lines and reproduces the resonance frequencies 
with a remarkable accuracy. The results of the calcula- 
tion are included in Fig. [2] and compared to the measure- 
ments in Tab.[l] The strongest molecular line, labeled by 
(0), in each spectrum corresponds to the dimer state lo- 
calized in the outermost potential well (Fig. [2]) which was 
assigned as vibrational ground state v — in [6 . The 
present study allows to identify several additional lines 
as excited dimer states. In particular, the theory pre- 
dicts a state which is even deeper bound than the state 
in the outermost potential well, (0)- This state is local- 
ized around R « 1400 ao for 35s (Fig. [5| and is assigned 
in the measured spectra. In addition to the dimer lines, 
we observe two resonances at deeper binding energy in 
the experiment. The lowest one systematically appears 
at twice the binding energy of the outermost localized 
dimer state (0). Our Green's function calculations show 
that for such large internuclear separations, the binding 
energy for two simultaneously bound ground state atoms 
is approximately given by the sum of the two binding 
energies. This approximation is a result of the weak in- 
teraction between the ground state atoms and the large 
geometric phase space of the trimer molecule. The ob- 
served spectral feature provides clear evidence for the di- 
rect photoassociation of a long-range triatomic molecule 
[T71 [18], in which two ground state atoms are simulta- 
neously bound to the Rydberg atom via their interac- 
tion with the Rydberg electron. The remaining unas- 
signed lines may indicate that excited states of these ex- 
otic trimers are also formed in our experiments. 

The agreement between theory and experiment con- 
firms existing s-wave scattering data, predicting a zero 
energy s-wave scattering length of A s (k = 0) = —16.05 ao 
[12J[T5]. In addition, the described fitting procedure al- 
lows to extract p-wave scattering parameters, and yields 
a zero energy scattering length of A p (k = 0) = — 21.15 ao, 
which is well within the range of previous calculations. 
Our extracted value for the shape resonance energy 
E res = 23.6 meV constitutes the first experimental deter- 
mination of this quantity and agrees well with theoretical 
predictions p~5j [19] . 



III. BINDING MECHANISMS 

For a pure s-wave potential the outermost localized 
state (0), for n = 35 at —23 MHz, constitutes the true 
vibrational ground state [6] of the molecule. A more ac- 
curate determination of the molecular potential requires 
the inclusion of p-wave scattering and in particular the 
shape resonance in the p-wave scattering cross-section of 
rubidium [I5j [19] has to be taken into account, which 
leads to an avoided crossing between adjacent molecular 
potential curves around ~ 1000 ao (see inset of Fig. 
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TABLE I: Observed and calculated molecular resonances. The frequencies are given in MHz. Errors (2a) of the measured 
binding energies are determined by the fit of the spectra. 



This has a dramatic effect on the vibrational structure 
and results in additional molecular states that are delo- 
calized over several 100 ao with vibrational energy that 
lies above the potential maxima inside the attractive po- 
tential well. Although energetically possible, these states 
have only negligible probability density at interatomic 
distances smaller than ~ lOOOao- As there exists no 
inner potential barrier that would facilitate a common 
molecular binding of these excited vibrational states, the 
binding has to arise from a fundamentally different mech- 
anism. We show that molecular binding in this case origi- 
nates from quantum reflection at the steep potential drop 
around ~ 1000 ao- 

Due to the enormous potential drop and flat potential 
energy curve at the avoided crossing, the molecular wave 
function behaves in the inner region very much like that 
of a free particle with a relatively large kinetic energy 
compared to the outermost part of the potential curve. 
The large discontinuity leads to a dramatic change in ki- 
netic energy over a small distance and the wavefunction 
is reflected at this point, similar to what happens at a 
large impedance mismatch in an electrical circuit, and 
the wavefunction piles up at large R. Consequently, the 
long-range molecular structure is largely unaffected by 
the details of the interaction at small distances, and inac- 
curacies in the shorter ranged energetically deeper parts 
of the interaction potential can safely be neglected. The 
small leakage into this region is treated as inward scatter- 
ing - in analogy to standard outward scattering problems 
- such that the molecular states can be obtained from a 
stabilization procedure [20]. Precisely, we place a reflect- 
ing boundary in the inner region of the potential well at 
Ro (see inset of Fig. [I]) and record the energy eigenvalues 
as a function of Ro. The observed energies are clearly 
visible in the resulting stabilization plot (see Fig(3^i) and 
appear as sharp maxima in the corresponding density of 
states (see Fig {3)3). 

Within an alternative approach that underlines the 
crucial role of quantum reflection, we consider an out- 
going plane wave injected at R = with an energy below 
the height of the innermost well. We study the complex 
reflection coefficient 1Z = \7Z\e lG of reflection off the outer 
potential plateau at R ~ 1000 ao- The scattering phase 



Q(E) can be related to the time delay r = dQ/dE be- 
tween a scattered and a freely propagating particle [21]. 
In Fig. [3jc) we show the phaseshift for the molecular 
curve of the n = 35 Rydberg state while its derivative 
is depicted by the right curve in Fig. |3|b). The energy 
dependence of r(E) exhibits a small background due to 
classical slow-down of the outgoing particle, moving up 
the potential hill. In perfect agreement with the calcu- 
lated density of states (left curve in Fig. [3^b)) we find 
an additional series of sharp resonances riding on top of 
this background. Consequently, these resonances are of 
non-classical origin and mark the occurrence of a quasi- 
bound state. The resonance line shape is of Breit-Wigner 
type, and, according to Wigner's time delay interpreta- 
tion, yields the lifetime of the corresponding quasi-bound 
state, in this case the time the system spends in the in- 
nermost part of the potential well before penetrating into 
the outer part of the potential. Adopting this to the vi- 
brational states in the outer part, r accords to the time 
the ground state atom spends in the outer region before 
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FIG. 3: (a) Stabilization diagram for n — 35, showing the 
vibrational energies as a function of the position of the in- 
ner, reflecting boundary at Ro (see FigQ. The corresponding 
density of states (dos), shown in (b), yields the energies and 
widths of the molecular resonances, in perfect agreement with 
the corresponding results obtained from the reflection phase 
(c). 
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entering the inner part of the potential. 

Such quasi-bound states, occur, e.g., in strong-field 
ionization of atoms and molecules [22 and affect trans- 
port phenomena in multi-layered semiconductors [23] , 
where localized states are confined inside a well by a high 
tunneling barrier. In contrast to these systems, the long- 
range vibrational states studied here arise nearly entirely 
from quantum reflection keeping the ground state atom 
outside the potential drop. To distinguish this purely 
quantum mechanics-based reflection from semiclassical 
phenomena, the spatial variation of the de Broglie wave- 
length A = 2it/^/2M[E - U(r)] has to be considered [24]. 
A a slow variation of A defines the semiclassical region. 
Depending on the vibrational energy, this condition is vi- 
olated, i.e. dX/dr > 1, around the potential maxima in 
the outer plateau region. Indeed the innermost of these 
'badlands' extends deep into the left of the sharp poten- 
tial drop. To identify the dominant contributions of the 
molecular curve, we have adapted the method of [25] to 
determine the reflection coefficient of the two inner po- 
tential maxima. We find, that even for the highest lying 
excited state for Rb(35s), the reflection coefficient is a as 
high as \K\ > 0.7. 

This clearly demonstrates that these unique vibra- 
tional states arise from the special form of the interac- 
tion potential leading to quantum reflection near the deep 
abyss due to the avoided potential crossing that results 
from the p-wave scattering resonance. 

IV. CONCLUSION 

The remarkably good agreement between experiment 
and theory over a range of principal quantum numbers 



attests to the accuracy of Fermi's original pseudopoten- 
tial approach in describing interactions in highly excited 
mult i- atom systems. A non-perturbative approach of 
these interactions is, however, necessary and leads to an 
unusual binding caused by internal quantum reflection. 
The latter is crucial for the existence of excited dimer 
states, since it prevents collapse to small interatomic dis- 
tances, which would lead to fast decay of the molecule. 
Since the molecular structure is found to be sensitive to 
the characteristics of both s-wave and p-wave scatter- 
ing, long-range Rydberg molecules are an accurate exper- 
imental platform to investigate electron-atom collision in 
a previously inaccessible ultralow-energy domain. 



The fact that these exotic molecules have huge bond 
length but nevertheless a vibrational spectrum with sev- 
eral states make a new type of ultracold Rydberg chem- 
istry feasible: larger polymers or even heteronuclear 
molecules can be realized, where the number of atoms, 
the constituent atomic elements and the addressed Ryd- 
berg state allow for the precise selection of the properties 
of the long-range molecule. Indeed, the three-body pho- 
toionization demonstrated in this work, may be regarded 
as an initiating step towards the production of larger and 
more complex polyatomic molecules. In particular, the 
high sensitivity of Rydberg atoms to external fields pro- 
vides an ideal starting point for coherent control and the 
investigation of molecular dynamics and reactions in a 
new regime. 
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V. METHODS 

A. Molecular potential calculations 

Our non-perturbative potential calculations are based 
on the short range character of the electron-ground state 
atom interaction. Hence, outside the range of V e A the 
electronic wavefunction 

i>(r) = J2CimG lm (r,E) (3) 

lm 

can be expanded in the solutions Gi m = 
{2i'-i)\\ [y^m(V / )G ? A(r, r', E)] r t=o of the atomic Hamilto- 
nian Ha- They are determined by the corresponding 
Green's function Ga, obtained with quantum defect the- 
ory. On the other hand, in the vicinity of the perturbing 
atom, we can neglect variations of the electron-core 
interaction [V^A+(| r — ~ V e js L +(R)]. In this case the 
Hamiltonian ([!]) describes a free-particle scattering from 
a finite-range potential [27], such that the wavefunction 

V(r) = ^c^l^n)^" 1 + ... + B^r 1 + ...)] (4) 

lm 

is entirely determined by the scattering phase shifts 
Si(k) via (21 - 1)!!(2Z + 1)!LB Z = k 2l+1 cot S t (k) of the 
different partial waves. Finally, the electron momen- 
tum k is expressed through the semi-classical relation 
E = k 2 /2+V e A+ (R), and its total energy E is obtained by 
matching the two different expansions Eq.([3| and Eq.Q 
as r — > 0. 



B. Experimental scheme 

The experiments are performed in an ultracold cloud 
(T = 3/iK) of magnetically trapped 87 Rb atoms in 
the ground state 5si/2,F — 2,rriF = 2. The nsi/ 2 
Rydberg state is addressed via a pulsed two-photon 
excitation: The first laser at 780 nm is blue-detuned 
with respect to the 5p 3 / 2 state and together with the 
second laser at 480 nm resonant to the transition 
5si/2 —> ns i/2- In the experiment, the two cw-lasers 
with a combined linewidth below 1 MHz are switched 
on for 3/is. For signal detection, we apply a field 
ionization pulse immediately after the laser excitation 
and probe the ions with a microchannel-plate. De- 
tails regarding the experimental setup and the sequence 
for the photoassociation spectroscopy can be found in [6] . 
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